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GTMDs in Light-Cone Quark Models 



Cedric Lorce 



1. Introduction 

Hadrons are composite objects. Their interaction with external probes like e.g. photons is pa- 
rameterized in terms of Lorentz-invariant functions. The most complete information is contained 
in the so-called generalized transverse-momentum dependent parton distributions (GTMDs). After 
appropriate Fourier transform, the GTMDs can be interpreted as Wigner or phase-space distribu- 
tions, giving access to the correlations between quark momentum and transverse position. 

A brief introduction to the generalized quark-quark correlator functions defining the GTMDs 
can be found in Section ||. On the light cone, this correlator can be written as the overlap of light- 
cone wave functions. We present in Section || the overlap restricted to the three quark (3Q) sector 
and specify the expression to a light-front constituent quark model (LFCQM) and the chiral quark- 
soliton model (#QSM). In Section |] it is explained why working on the light cone is mandatory 
to develop a (quasi-)probabilistic interpretation of the distributions. Finally, we apply in Section || 
the formalism presented in the previous sections to study the distribution of an unpolarized quark 
in an unpolarized proton. 

2. General Quark- Quark Correlator 

The maximum amount of information on the partonic structure of the nucleon is contained in 
the fully-unintegrated quark-quark correlator W for a spin- 1/2 hadron [jj], ^, ||, f|], defined as 

W^ A (P,k,A,N;r 1 ) = y t|^ z (p' ,A'\W(- L 2 z)rW ¥ (±z)\p,A). (2.1) 

This correlator is a function of the initial and final hadron light-cone helicities A and A', the average 
hadron and quark four-momenta P = (p' + p)/2 and k, and the four-momentum transfer to the 
hadron A = p' — p (see Fig. |l]for the kinematics). The superscript T stands for any element of the 




P-A/2 P + A/2 



Figure 1: Kinematics for the fully-unintegrated quark-quark correlator. 

basis {1,75,7^, 7^75, Z(7^ v } in Dirac space. A Wilson line W = W(—^z,^z\n) ensures the color 
gauge invariance of the correlator, connecting the points — ^z and ^z via the intermediary points 
— \z + 00 • n and ^z + 00 • n by straight lines. This induces a dependence of the Wilson line on the 
light-cone direction n. Since any rescaled four-vector an with some positive parameter a could be 
used to specify the Wilson line, the correlator actually only depends on the four- vector N = y^-, 
where M is the hadron mass. The parameter X] = sign(n°) gives the sign of the zeroth component 
of n, i.e. indicates whether the Wilson line is future-pointing (17 = +1) or past-pointing (tj = — 1). 
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The quark-quark correlators parametrized in terms of generalized parton distributions (GPDs), 
tranverse-momentum dependent parton distributions (TMDs) and form factors (FFs) correspond 



to specific limits or projections of Eq. (|2JJ). These correlators have in common the fact that the 
quark fields are taken at the same light-cone time z + = 0. Let us then focus our attention on the 



k -integrated version of Eq. (2.1) 



W$ A (P,x,k±,A,N;ri) = J dk-W^ A (P,k,A,N;r]) 



1 / ^^^ z {p'^M-kz)TW W (\z)\pA) 



(2n 



(2.2) 



z +=o' 



where we used for a generic four- vector = [a + ,a~ ,a±] the light-cone components a ± = (a ± 
a 3 )/\^2 and the transverse components a± = (a 1 , a 2 ), and x = k + /P + is the average fraction of 
longitudinal momentum carried by the quark. A complete parametrization of this object in terms 
of so-called generalized transverse momentum dependent parton distributions (GTMDs) has been 
achieved in [Q]. The GTMDs can be considered as the mother distributions of GPDs and TMDs. 
For example, the (T-even) distribution of unpolarized quarks in an unpolarized hadron is given by 
the GTMD F( x , and is related to the GPD H and the TMD f x as follows 

tf(x,0,Ai)= / d\ L F e n {x^kl,k L -A L: A\), 

J - (2.3) 

f [ (x,kl)=F^ l (x,0,ki,0,0). 



3. Overlap Representation 

Following the lines of []|, g], we obtain in the light-cone gauge A + = an overlap representa- 
tion for the correlator (2.2) at the twist-two level restricted to the 3Q Fock sector 1 



1 



wP A (P,x,k±,A,N;r 1 )= 



£ / [dx] 3 [d 2 fc J 3 A{k) Yxp (/) Yap (r) U M ^ . 



(3.1) 



i=i 



where the integration measures are defined as 



[dx] 3 



3 

JJdXi 

i=\ 



5 , [d 2 ^] 3 



f=i 



3 d 2 k 



n 2 (27r)3 



KlnfS® [fk iX . (3.2) 



v i=l 



Furthermore, in Eq. ( ^1| ) the function A(k) = 30(xi) 5(x— xi) d^(kx — ki±) selects the active 
quark average momentum (we choose to label the active quark with i = 1 and the spectator quarks 
with i = j = 2,3). The 3Q LCWF YAj}( r ) depends on the momentum coordinates ki = (y,-, fQjJ of 
the quarks relative to the hadron momentum (collectively indicated by r), and the index j8 which 
stands for the set of the quark light-cone helicities {A,-}. The transition from the initial quark 
light-cone helicity A, to the final one X( is described by a complex- valued 2x2 matrix M x>i ^. 
In particular, we have for the spectator quarks M j > = o ■> '. For the active quark, the matrix 



'Quark flavor and color indices have been omitted for clarity. In the processes considered here the flavor and color 
of a given quark remain unchanged. 
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^fAj'Ai depends on the twist-two Dirac structure r tw i st _2 = {y + ,y + 75,JCT +1 y5,/a +2 y5} use( l in the 
correlator, see e.g. [0, g]. 

We choose to work in an infinite momentum frame such that P + is large, P± = 0± and A • P = 0. 
The four-momenta involved are then 



P 

k 



p+ m t 4 



'2(1-£ 2 )P+' 

xP + ,k~,k ± 



-25P+,5 



0,±1,0± 



M 2 + ^ x 
;i-^ )jP+ > - 



(3.3) 



Note that the form used for « is not the most general one, but leads to an appropriate definition 
of TMDs for SIDIS and DY processes. For the active and spectator quarks the initial and final 
momentum coordinates are then 



x + 5 1 -x Aj_ 



1 + 5 



Xi A I 



3 



1-5 



1-5 



Y/'J- 



x; A ± 



(3.4) 



1-5 2 



So far, the exact 3Q LCWF derived directly from the QCD Lagrangian is not known. Nev- 
ertheless, we can try to reproduce the gross features of hadron structure at low scales using con- 
stituent quark models. Many models exist on the market based on the concept of constituent quarks. 
However only a few incorporate consistently relativistic effects. We focus here on two such mod- 
els: the light-front constituent quark model (LFCQM) [[7], ||, ^] and the chiral quark-soliton model 
(#QSM) ijfij |l|, II, HI, |4|. The LCWFs used in LFCQM and in #QSM have a very similar 
structure given by 

¥A P (r) =^(r)£or 2ff3 lW.(% (3-5) 

Oj i=l 

where *P(r) is a global symmetric momentum wave function, <J>^ lCT2f73 is the SU (6) spin-flavor wave 
function, and Dik) is an SU (2) matrix connecting light-cone helicity A,- and canonical spin a, 



D(k) 



1 

\K\ 



~ K R K z , 



K R ,=K l ±iK 2 . 



(3.6) 



The explicit expressions for the momentum wave function *P(r) in Eq. (3.5) and the vector K in 
Eq. (0) in LFCQM read 



¥(#■) =2(27i) 



(O1CD2CD3 



xix 2 X3-#o (X) 2 + i8 2 ) r ' 



X, = m + y^ , Kx = Kx, K z =y^ — (O, 

(3.7) 

where jV is a normalization factor, ^#0 = L; ^ i s tne f ree invariant mass, ft),- is the free energy 
of quark i, m is the constituent quark mass, and j8, y are model parameters fitted to reproduce the 
anomalous magnetic moments of the nucleon [|l5|]. On the other hand, within the #QSM one has 



y{r)=jrY[\Ki 

(=1 



Kz = h + j4rj, 



K 7 



(3.8) 
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where is the soliton mass, E\ ev is the energy of the discrete level in the spectrum, and h, j are 
the upper and lower components of the Dirac spinor describing this discrete level. 
For further convenience we introduce the tensor correlator 



Z A'A 



(3.9) 



where W£, A = (wj$, W^ >75 \ wJ^^W^ 1 ) and 6^ = (l,d) with a ; the Pauli matrices. We 



now use the LCWF given by Eq. (3.5) and write the overlap representation of the correlator tensor 
W» v as 

W» v (P,x,kt,A,N;r,) = J_ / [dx] 3 [d 2 *j_] 3 A(k)W*(r')^(r)^ v (r' ,r), (3.10) 

a/1 —c 1 J 



where stf^i/ \r) stands for 



v (r',r) =AO$ v (h -h)+B [l$ (Z 3 • 0{? + 1% (Z 2 • Oi) v ] • 



(3.11) 



In Eq. ( |3.11| ), Zy = and the matrix is given by 



^y 



1 



l^'ll^l 



i K'xK 



i K'xK 



i K'xK 



K'K 

i (K'xK) K'K- 2K' x K x —K' x K y - K' y K x K' X K Z + K[K X 
i (k' x K\ —K' y K x - K' x K y K'K- 2K' y K y K' y K z + K' 7 K y 



V 



i ( K' x K J -JSjK, - JSjiQ - -K' K + 2K[K 



(3.12) 



The tensor correlator W^ v in Eq. (3.9) has two indices. The index /I refers to the transition in terms 
of baryon light-cone helicity, while the index v refers to the transition in terms of the active quark 
light-cone helicity. For example, the components W and W 03 correspond to the matrix elements 
of the 7 + and 7+75 operators in the case of an unpolarized hadron, respectively. Equation fl3.1C ) 
gives the explicit expression for the tensor correlator in terms of the overlap of initial 'P(r) and 
final *F*(r') symmetric (instant-form) momentum wave functions with the tensor £/ tlv (r J ,r) for 
a fixed mean momentum of the active quark k. The tensor £/^ iV (r J ,r) contains the spin-flavor 
structure derived from the overlap of the three initial and final quarks. Taking into account the 
possible couplings of the helicities of the active and spectator quarks to give the hadron helicity, 



the coefficient A and B in Eq. ( p.l 1| ) for SU(6) spin-flavor wave functions are 



A p 



1, 



2. 



(3.13) 



Furthermore, the matrix 0^ v in Eq. (3.12) describes the overlap of the initial and final quark states. 
The columns are labeled by the index v which indicates the type of transition in terms of quark 
light-cone helicity. The rows are labeled by the index /I which indicates the type of transition in 
terms of quark canonical spin. This matrix reduces to if = Of for the spectator quarks, since in 
this case the light-cone helicity is conserved. 
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4. Impact Parameter Space 



According to the standard interpretation [17, 18], the charge density can be identified with the 



three-dimensional Fourier transform of the electric Sachs Form Factor Ge, i-e. the matrix element 
of the current density 7° 



This identification is actually only valid in the nonrelativistic approximation. To work out the 
Fourier transform, one has to know the form factors for every Q 2 . In the Breit frame the latter is 
identified with the three-momentum of the virtual photon Q 2 = q 2 . This means that for every value 
of Q 2 we have to move to a different frame and the charge density undergoes naturally a different 
Lorentz contraction. Moreover, in order to have a probabilistic/charge density interpretation, the 
number of particles should be conserved. However, in the Breit frame nothing prevents the photon 
to create or annihilate a quark-antiquark pair. 

All these problems are cured in the infinite momentum frame with q + = (the so-called Drell- 
Yan-West frame). The photon is kinematically not allowed to change the number of quarks as 
the light-cone momentum of a massive particle is strictly positive p + > 0. Moreover, the hadron 
undergoes an extreme Lorentz contraction and looks like a pancake. Only a two-dimensional charge 
density [jjj, is then meaningful and can be identified with the two-dimensional 



Fourier transform of the matrix element of J 



+ 



dV 

(1%) 2 2p+ 



P3<P±)= I 7^7 -^(p + ,^,s\J + (0)\p + -^,s), (4.2) 



the photon virtuality being given by Q 2 = q 2 , and s denoting the hadron polarization. 



5. Wigner Distributions 

Wigner distributions are quantum phase-space distributions, containing all the correlations be- 
tween position and momentum of the partons. Since Heisenberg's uncertainty principle forbids to 
determine precisely both position and momentum of a quantum state, Wigner distributions have 
to be considered as quasi-probabilistic distributions. In the context of quantum field theory, they 
have already been discussed to some extent in the Breit frame [pp. For the reason mentioned in 
the previous section, it is actually preferable to work in the infinite momentum frame. By Fourier 
transforming in Eq.(^3|) with respect to Aj_ and integrating over x, we obtain a (transverse) 
phase-space distribution p(k 2 L ,k± ■b± 1 b 2 L ). The only two available transverse vectors are k± and 
b±. This means that for fixed kx • b± and \k±\ the distribution is axially symmetric, which is physi- 
cally meaningful since there is not a preferred direction in the transverse plane and a global rotation 
around e z should not have any effect on the distribution of unpolarized quarks in an unpolarized 
nucleon. 

If we are now interested in the amplitude to find a fixed transverse momentum k± in the trans- 
verse plane, the distribution is not axially symmetric anymore due to the k± • b± dependence. Al- 
though the Wigner distributions can not be directly extracted from experiments, they can be inferred 
from the LCWFs parametrized on the basis of our phenomenological knowledge of GPDs, TMDs 
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Figure 2: Typical distribution (left) and equi-amplitude line (right) of fixed k ± — k ± e y in the transverse 
plane for an unpolarized quark in an unpolarized nucleon. 



p [l/GeV 2 fm 2 ] 



1/5, 


% \ 


/ * 


* \ 


/* 


*v 




& 


J 0,5 





-1.0 -0.5 



0.5 1.0 



b [fm] 



Figure 3: Typical radial distributions in the transverse plane for fixed k± ■ b and \k± | of an unpolarized quark 
in an unpolarized proton. Solid line corresponds to half of the up quark distribution. Dashed line corresponds 
to the down quark distribution. 

and form factors. As a starting point, we used the ^QSM and LFCQM which were introduced in 
Section || and already tested in the calculation of several nucleon observables [12, [l3|, 14, 15, |l6| ]. 



In the case of Ffj, both models give the same qualitative picture with a larger distribution 
amplitude when k± _L b± and smaller when k± \\ b± (see Fig. This can be understood with 
naive semi-classical arguments. The radial momentum (k± -b±)b± of a quark has to decrease 
rapidly in the periphery because of confinement. The polar momentum k± — (kx -b±)b± receives a 
contribution from the orbital motion of the quark which can still be significant in the periphery (in 
an orbital motion, one does not need to reduce the momentum to avoid a quark escape). This naive 
picture also tells us that this phenomenon should become more pronounced as we go to peripheral 
regions (\b±\ 3>) and to high quark momenta (\k±\ S>). This tendency is supported by both models. 

It is also interesting to compare up and down quark distributions. Let us fix once again k± ■ b± 
and |jfcj_|. The distribution being axially symmetric in the transverse plane, we focus on the radial 
distribution, see Fig. |3| The up quark distribution has been divided by two for comparison with 
the down quark distribution. Up quarks appear to be more concentrated around the center than 
down quarks. For a neutron, we just have to exchange up and down quarks. We can therefore 
see that the center of the neutron is negative, in agreement with the conclusion obtained using the 



phenomenological neutron FFs [22] 
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6. Summary 

We used the framework of light-cone wave function to study generalized transverse-momentum 
dependent parton distributions which parametrize the most general quark-quark correlator. These 
distributions are connected to the Wigner or phase-space distributions by Fourier transform. We 
presented a general expression for the 3Q contribution to the generalized correlation functions 
defining the GTMDs and applied it to a light-front constituent quark model and the chiral quark- 
soliton model. A non-trivial pattern for the phase-space distribution in the transverse plane for 
an unpolarized quark in an unpolarized proton has been observed and interpreted semi-classically 
as related to the quark orbital angular momentum. We also confirmed the picture of a negatively 
charged core in the neutron. A presentation of the systematic study of generalized transverse- 
momentum dependent parton distributions is in preparation. 
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